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Abstract 



D ■ 

Geometrical CP violation is a particular type of spontaneous CP violation in which 
the vacuum expectation values have phases which are calculable, i.e. stable against the 
£N| \ variation of the free parameters of the scalar potential. Although originally suggested 

within a specific version of the three-Higgs-doublet model, it is a generic phenomenon. 
We investigate its viability and characteristic features in models with several Higgs dou- 
\Q \ blets. Our work contains both general results and illustrative examples. 

m 

^ ; 1 Introduction 

1.1 Geometrical CP violation 

CP violation [lj is a hallmark feature of the electroweak interactions. Despite a long his- 
tory of experimental and theoretical exploration, its fundamental origin remains enigmatic. 
Various models of physics beyond the Standard Model (bSM) try to provide an explanation 
of that origin. One attractive idea is that CP is a valid symmetry of the Lagrangian but 
gets broken spontaneously by the vacuum state. In a seminal work, T. D. Lee [2] has shown 
that spontaneous CP violation can take place in the scalar sector of the two-Higgs-doublet 
model (2HDM); later on, when the existence of three quark generations became known, it 
was shown |3j that the phase generated by spontaneous CP violation in the vacuum state of 
the 2HDM actually propagates into the quark mixing matrix. Since then, spontaneous CP 
violation in the 2HDM remained a rich subject, with several of its facets uncovered just quite 
recently [I]; a useful review is ref. [5]. 

In the 2HDM, spontaneous CP violation arises in the form of a relative phase between 
the vacuum expectation values (VEVs) of the neutral components of the two Higgs doublets, 
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despite the fact that all the parameters of the scalar potential are realJ3 This relative phase de- 
pends on the numerical values of those parameters. This fact renders the model less predictive 
than one might have wished, since one can tune the parameters of the potential and thereby 
produce any phase that one wants; in addition, the phase is subject to renormalization. 

With larger scalar sectors, a new form of spontaneous CP violation known as geometrical 
CP violation can take place. In it, the relative phases of the VEVs become calculable, i.e. 
their values do not explicitly depend on the free parameters of the potential, rather they arise 
from its structural properties!! Besides having a more transparent origin, such phases are also 
stable against renormalization. 

Geometrical CP violation was originally discovered [7j in the three-Higgs-doublet model 
(3HDM). Generically, the tree-level scalar potential of that model can be written as a sum of 
one part which is invariant under any phase rotations of the doublets, 



N r 21 
v = J2 - m l<t>i<Pk + h [4>l4>k) + Yl [ Xkl {^<Pk) (frfi) + K 



k=l 



Kk<KN 



with N = 3, and a phase-dependent part V^hase- I n the model considered in ref. [7J, the poten- 
tial was assumed to be invariant under the finite group A (27) of Higgs-family transformations 
which is generated by the cyclic permutations of the three doublets and by the rephasing 
transformation 

02 W0 2 , 03 ->■ W 2 3 , (2) 

where u = exp (i 2iv/3). With this A(27) symmetry, the scalar potential is very restricted and 
only one phase-sensitive term survives: 



^phase A 



^02 ) 0l03 + 0203 0201 + ( W ( W + H.C. . (3) 



It was noted in ref. [7J that, for positive A, the alignment of the VEVs at the global minimum 
is either 

((ft.ift-fft)^^^), (4) 

with positive real u, or other alignments related to the one in eq. (jl]) by complex conjugation 
and/or by a transformation belonging to the symmetry group A (27). 

The subject of geometrical CP violation lay dormant for decades but, recently, interest 
in it has revived. Continuing the analysis of the original A(27)-symmetric 3HDM potential, 
the authors of ref. [8J have shown that phases remain calculable even when effective non- 
renormalizable scalar interactions of arbitrarily high order are present in that model. In ref. [H] 
the idea of geometrical CP violation was extended to models with several Higgs doublets. On 
a more mathematical side, a recent work [10] has argued that calculable phases, and in general 
any stable kind of alignment of VEVs, may arise in models whose orbit space has sharp corners; 
this remark provides an additional perspective on the meaning of the word 'geometrical' in 



1 A different form of spontaneous CP violation, in which it occurs in spite of all the VEVs being real, has 
been proposed in ref. [6]. 

2 The statement that the relative phases of the VEVs in geometrical CP violation do not depend on the free 
parameters of the potential is sometimes found. This statement is somewhat misleading since those parameters 
must be such that the CP- violating VEV alignment is the global minimum of the potential. A better definition 
of calculable phases is that they are phases which are stable against reasonably large variations of the free 
parameters of the potential. 
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the phrase 'geometrical CP violation'. Finally, in ref. [TT] the question of when calculable 
phases lead to genuine CP violation was recast in the form of a distinction between inner 
and outer automorphisms of the Higgs-family-symmetry group of the potential; that work, 
however, referred only to the original 3HDM example of geometrical CP violation in ref. [7]. 



1.2 Goals and scope of this work 

In this work, we report on the results of a study of geometrical CP violation in multi-Higgs- 
doublet models that we have undertaken. 

Our work differs from ref. [9] in several aspects. Firstly, after deriving calculable phases 
of VEVs, one has to check that they cannot be eliminated by a rephasing symmetry of the 
potential. This check becomes automatic once we know the rephasing symmetry group of 
the given potential. So, we start our investigation by developing a theory for finding it for 
any iV-Higgs-doublet model (NHDM) potential. Secondly, we consider only renormalizable 
potentials, since they offer already quite a rich spectrum of possibilities. The inclusion of 
non-renormalizable terms goes beyond the scope of this work, but they can also be analyzed 
using the methods developed here. Thirdly, we limit the variety of cases by requiring that the 
phase-sensitive part of the potential contains only one cyclic chain of terms, which is protected 
by a specific rephasing symmetry group not shared with any other possible terms. We call 
such a symmetry 'stable' because it guarantees that no other renormalizable terms in the 
potential are generated upon renormalization. Finding out which rephasing symmetry groups 
are stable and what are the potentials that they stabilize is another task that we solve in this 
paper. In section H] we shall also comment on the consequences of lifting this requirement. 

In order to ensure that no additional terms arise in the potential upon renormalization, 
one usually imposes a symmetry on the Lagrangian. In the particular case of geometrical CP 
violation, it is known [7] that one must impose an additional symmetry beyond just rephasing, 
because otherwise one would get only trivial calculable phases or 7r among the VEVs. In the 
original work [7], this was formulated as a requirement of a non-Abelian symmetry group of 
the potential, but this formulation is ambiguous because not every Abelian symmetry group 
of a NHDM potential can be represented by rephasing transformations [13J. 

Specifically, we consider a NHDM whose Higgs potential is invariant under the Z^v group 
of cyclic permutations of the doublets 

01 -> 02 -> > <j>N -> 01, (5) 

and also under the CP transformation <p k (t,x) — > <fi* k (t,—x) for k = 1,2, ...,N. These 
invariances significantly restrict the form of Vo but still leave much freedom in Vphase- Indeed, 
Vphase may contain both quadratic terms of the form 

V k = \k [4>\(f>k + 020fc+i H 1- 0jv0fc-i + H.c.J (6) 

and quartic terms of the form 



AM (0?0m) + (020fc+l) (0l + i0m+l 



+ •■■+( 0t0fe-l ] ( 0Ll 0m-l ) + H.C. 



(7) 
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In eqs. flS]) and ([7]), and below, the subscripts are always understood to be modulo N. Since 
the Z N cyclic symmetry (jSJ) will always be present in our examples for Vphase, we shall from 
now on, for the sake of brevity, write down only the first term followed by the expression 
"+ cyclic." 

The analysis of spontaneous CP violation in the NHDM with a generic collection of terms 
Vk and Vkim is a complicated task which goes much beyond the scope of this paper. Here, we 
just want to present some simple examples realizing the idea of geometrical CP violation for 
N > 3. For this purpose, we impose an additional invariance of the potential V^hase under a 
certain rephasing symmetry group G; that additional invariance selects only one set of phase- 
sensitive terms out of all the Vk and Vkim- We call the rephasing symmetry G stable if it 
guarantees that Vphasc contains only one cyclic set of renormalizable terms compatible with it, 
so that no additional symmetric terms may appear upon renormalization. 

The criterion of stability immediately allows us to disconsider all the terms 14 in eq. flS}. 
Indeed, if V& were invariant under a certain G, then Vkik would also be invariant, and this 
disagrees with the definition of stable symmetry. Thus, in this paper we only have to deal 
with quartic terms V&j m . 

Finding out stable rephasing symmetries is our first task and is solved in section [2j When 
one such symmetry is found, we then need to minimize the potential and check whether 
geometrical CP violation may indeed take place in the global minimum. 

Minimization of a NHDM potential with several terms is a highly non-trivial problem and 
we do not attempt to solve it in full generality. We merely argue that it is always possible to 
find a region in the space of parameters such that the VEVs are of the form of pure phases, 
viz. (0°) o = wexp (i6k) with some positive real u which is the same for all k. We then need 
to find the values of the 9k at the global minimum. If those phases are non-trivial and if, 
moreover, they cannot be all rotated away through a symmetry of the potential, then we 
have spontaneous CP violation. In order to be able to call that CP violation geometrical, 
the phases 9k must be calculable, i.e. they cannot explicitly depend on the parameters of the 
potential. 

2 Stable symmetries 

2.1 The rephasing symmetry group and the Smith normal form 

The first question that one faces is how to determine the rephasing symmetry group of a given 
Vkim- This question can be answered in its full generality by using the methods developed in 
ref. [13]. (See also refs. [12] , where a similar construction has been used.) Here we slightly 
adapt those methods to our problem. Some basic elements of the machinery that we are going 
to expose in this section were also used in ref. [9]. 

Let us consider the full rephasing symmetry group [{7(1)]^, in which the k-th U(l) factor 
stands for phase rotations of k-th Higgs doublet, (pk — > 0fcexp (ia&). Consider one particular 
quartic term in Vki m , for example the first term (0[0fc)(0^0 m ). Upon the action of [{7(1)]-^, 
this term acquires a phase shift 



N 
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This term is invariant if £1 = 27rni with some integer n\. Similar conditions, viz. = 2irnk 
with integer n& for all k < N , must hold for all the quartic terms in Vki m - We thus obtain the 
system of equations 

N 

^ d kj aj = 2irn k , (9) 
j'=i 

which is expressed via anJVxiV matrix of integers -D = [dfcj]. The entries of D are either 0, 
±1, or ±2, and its rows may only be — up to the overall sign and to cyclic permutations — of 
one of the three types 



„•••>*'> 
1,... 



1 1 

2,... 



2 —2 
1 



-1,...") 



-l,...,±l,...,Tl,. 
where the dots indicate zeros. For example, for 



Vo 



213 



we have 



A 

/ 2 




218 



Ah 



103 ) + cyclic + H.c. 



(10a) 
(10b) 
(10c) 



(11) 



D 



\ 








1 


















2 -1 
2 




\ 


-1 
-1 

2 J 



(12) 



i.e., the rows are of type ( llObj) . 

The solutions of eqs. ([9]) form the rephasing symmetry group G of V^im- There is a sys- 
tematic procedure for finding out this group [15] . It is known that by performing elementary 
transformations of the following types: 

adding one column or one row to another column or row of D, 

flipping the sign of one column or one row of D, 

permuting columns or rows of D, 

one can bring D to its Smith normal form, which is a diagonal matrix 

diag (g?i, d 2 . . ■ , d r , 0, . . . , 0) (13) 

of positive integers c4 such that each dk is a divisor of dk+i- In eq. (fT3|) . r = rank D. Since the 
elementary transformations described above preserve the structure of eqs. we thus obtain 
a system of uncoupled equations 

271-%, (14) 



for j = 1, 



r, where the a,- are linear combinations of the initial a^, and the fL- are integers. 



The rephasing symmetry group of Vkim is then 



G = Z dl x Z d2 x • ■ ■ x Z dr x [U(l)] N - r 



(15) 
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where Zi is the trivial group. This result solves in principle the problem of finding out the 
rephasing symmetry group of any Vki m - 

Some information on G can be extracted even without explicitly computing the Smith 
normal form of D. Since all the elementary steps used in the computation of the Smith 
normal form preserve the absolute value of det D, one concludes that either |det D\ = Ylk=i ^ k 
if r = N or det D = if r < N. In the case of the scalar potential of a NHDM, Vki m is 
always invariant under one and only one U(l) group, which is the gauge group of overall 
phase rotations (pk — > 0fcexp(ia); no other U(l) symmetry group may exist lest Goldstone 
bosons arise upon spontaneous symmetry breaking. We therefore conclude that the matrix 
D is always singular with r = N — 1. The invariance under the U(l) gauge group does 
not represent any interesting property of the potential. We may remove this uninteresting 
rephasing symmetry by considering the (N — 1) x (N — 1) matrix D which is obtained from D 
by deleting a single column and a single row of D. (These may be any row and any column.) 
The remaining part of the rephasing symmetry group, G/U(l), is finite and has order 

\G/U(1)\ = detD . (16) 

Now, if, in particular, it happens that | det D\ is a number whose prime decomposition involves 
only primes with power one, i.e. \ detD\ = pip 2 ■ • -p s , then the rephasing symmetry group is 
uniquely determined to be 

GMl)=Z |dct5| , (17) 

since in that case Z Pl x Z P2 x • • • x 7L Vs is the same as Z| dot ^,|. 

In other cases, this simple procedure is not sufficient to find the rephasing symmetry group 
and one needs to actually compute the Smith normal form. Finding the Smith normal form is 
easily algorithmizable and there exist computer-algebra packages [H] which do it for any input 
matrix of integers; these packages are adequate for practical purposes and for any reasonable 
number of doublets. 

To keep the notation short, we shall from now on suppress the U(l) factor in the rephasing 
symmetry group G, implicitly assuming that it is present everywhere. 

2.2 The rephasing symmetry group of V213 with generic TV 

Consider the prototypical potential V213 in eq. ( ITT]) for a generic number of doublets N. Let 
us suppose that G/U(l) keeps the first doublet <\>\ invariant and that it transforms the second 
doublet 02 with a phase factor a. We furthermore suppose that the k-th doublet is transformed 
as <pk o" Sfc 0fc, with s\ = and s 2 = 1 by construction. Using the form of the potential, 
which contains terms (0fc_ 2 0fc-i)(0k_20fc)j one can derive a recurrence relation for the Sk, 

s k = 2s fc _ 2 - (18) 

Solving this recurrence relation via the Ansatz s k = a k , we obtain that a may be either 1 or 
—2. Inputting the initial conditions s± — and s 2 = 1, one finally obtains 

1 - (-2) h - 1 , s 

Sk = , (19) 

which gives the sequence 0, 1, —1, 3, —5, 11, —21, 43, —85, 171, . . .. 
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The fact that the rephasing transformation is generated by a single factor a means that 
the group G is in this case of the form Z ? . In order to find the value of q, we note that the 
doublet (pN+i is identical with (pi, and therefore must transform with the same phase factor 
1, viz. <j Sn+1 = 1. Therefore the rephasing symmetry group of V213 is 

G 213 = Z\ SN+ll . (20) 



2.3 Stable symmetries 

The next question that we ask is whether Vum has a stable rephasing symmetry group, i.e. a 
rephasing symmetry group that allows only Vki m in Vphase- 

We firstly note that Vki m and V m ik transform in the same way under a rephasing transfor- 
mation. Therefore, they have the same rephasing symmetry group. Since a stable symmetry 
must allow only one cyclic set of terms, we conclude that Vki m must be identical with V m ik- 
This is possible only when I = 1 or (which is equivalent) when k = m. As a consequence of 
this criterion, from now on we shall consider only sets of terms of the form Vki m , i.e. with 
/ = 1. 

Next, consider the matrix D for some Vk\ m - If it turns out that by summing or subtracting 
rows of D one can obtain another row — distinct from all the rows of D — of one of the types (flOj) . 
then the symmetry is not stable. If this is impossible, then G is stable. 

We illustrate this criterion by considering V213 for N = 4. We know from the previous 
subsection that then G = Z 5 . The matrix D is 



D 



/ 2 -1-1 \ 
2-1-1 

-10 2-1 

\ -1 -1 2 J 



(21) 



The sum of the first and third rows of this matrix gives (1, —1, 1, —1), which is of the 
type fllOcp and, thus, allowed. This means that the terms (0^0 2 )(0i03) and (0304)(030i), 
which are present in V213 for N = 4, generate another quartic term (0i 02) (0304), which is also 
invariant under the same Z5 but is absent in V213. Therefore, V213 is not protected by a stable 
G for N = 4. In fact, it can be shown that in the case N = 4 there is no stable G for any 

^phase • 

In the case N = 5, on the other hand, the sum of the first and third rows of D yields 
(2, —1, 1, —1, —1), which is not an allowed row. In this case the terms (^j^X^I^s) and 
(0304) (0305) generate only non-renormalizable effective terms with six doublets, but no other 
renormalizable terms. One may check that no other combination of rows leads in this case 
to a new allowed row. This means that for iV = 5 the potential V213 is protected by a stable 
symmetry, which from the previous subsection we know to be Z n . 



2.4 Stable symmetries for small iV 

We apply here the methods described above to list all the potentials Vki m invariant under 
stable symmetries for small numbers of doublets. 

For iV = 3 there is only V213, which is invariant under a stable symmetry Z 3 and was 
studied in the seminal paper of ref. [7]. 
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For N = 4, as we said above, there is no Vum which can be protected by a stable symmetry 
group. 

For N = 5, only 

V213 , V315 , V412 , V514 (22) 

are protected by a stable symmetry, which is Zn for all the potentials in (1221) but with different 
generators. Generically, the potentials in ( 1221 can be written as Vi+k, 1,1+2* with fc = 1,2,3, 4 
(we remind the reader that all subscripts are meant modulo 5 in this case). In each of these 
four cases, the generators acts on the doublets in the following way: 

01 ->■ 01, 01+fc -> O-01+fe, 01+2A; -> 0" 1 Vl+2fc, 01+3fc O" 3 01+3fe, 01+4fc -> O" 6 01+4fc, (23) 

where a 11 = 1. 

We do not consider the case N = 6 in detail here because, as we shall see later, there is no 
geometrical CP violation for stable symmetries in the case of even N. 

For N = 7 there are two distinct sets of potentials, with six Vum in each set. The generic 
form of each set and its respective symmetry group are 

Vi+M,i + 2k with G = Z 43 and K+k,i,i+3fc with G = Z 29 , (24) 

for k = 1,2,3,4,5,6. Again, within each set of V^i m , the generators realizing the symmetry 
group are different for different potentials. For the second set in (T241 . the symmetry group G 
is 

01 ->■ 01, 01+fc -> P01+fc, 01+2fc -> P 8 01+2fe, 01+3fc -> /O 28 01+3fc, 

01+4A: - > p </>l+4fc, (Pl+Sfc — > P <Pl+5k, <Pl+6k ~^ P <Pl+6k, 

where p 29 = 1. 



3 Calculating the phases 
3.1 Minimizing the potential 

The simplest form of geometrical CP violation is characterized by VEVs which differ only by 
phases: (0^) o = ue l6k . Before calculating the phases 6k, we need to prove that the potentials of 
the type that we consider indeed allow, without finetuning, for a global minimum with equal 
|(0°) o | 2 = u 2 for all k = 1,...,N. 

Let us first consider only the potential Vq. We assume that its the minimum is neutral, i.e. 
that it does not violate electric-charge conservation. Substituting the VEVs in Vq, we rewrite 
it in terms of u\ = |(0fc) o | 2 - The expression must respect the cyclic symmetry and therefore 
it must be of the form 

N ^ N 

(Vo) = -«oE^ + 2E A ^ M < • ( 26 ) 

k=l k,l=l 

The real and positive-definite matrix is symmetric and circulant: 

A w = A ik and A fe +i,i+i = A H , Vk,l = l,...,N. (27) 
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Each row of Ay is a cyclic permutation of the previous row, and therefore the sum of all 
elements in a single row, A = Ay, is the same for all rows. Equation ff26l) can then be 
rewritten as 

W„4E A «K-4K-<)-^, where U g = ^. (28) 
lfe,i=i 

It then becomes evident that the minimum of Vq occurs, when m§ is positive, at u\ = Uq, viz. 
the minimum respects the cyclic symmetry Zjv- 

Next we need to prove that the inclusion of the phase-sensitive terms in Vki m does not spoil 
this cyclic symmetry and only amounts to a shift Mo — >* u - Upon substitution of the VEVs 
(02)o = u fc e ** fc > ^fcim becomes 

TV 

(Vfcim)o = 2 -^fcim (u\u k Um cos ^i + cyclic) , where ip k = ^ d w 0j, (29) 

z=i 

with the same matrix D as in section 12.11 In order that uq evolves smoothly to u while keeping 
all the Uj equal, we need to require the partial derivatives d (14i m ) /duj to be equal for all j. 
Computing those partial derivatives, we observe that we need to require that the quantities 
Cj = 2 cos ipj + cos ipj + k~i + cos ^j +m _i are all equal at the minimum. This happens only when 
all the cosipj are equal. Later we will see that this condition is indeed verified at the global 
minimum, which confirms that the Ansatz (0°) o = ue tBk is valid. 



3.2 Calculating the phases: general theory 

In this section we consider Vkim and ask what are the phases of the VEVs, viz. the 9k, that 
will minimize it. At the minimum, the expectation value of V^xm is 2Afci m u 4 J, where 

J = COS^i + COS ^2 + • • ■ + COSIpN- (30) 

Obviously, the stability points of J are points where 

dJ N N 

— = - sin i>i d i>t = ~Y1 (° T ) kj sin ^ = ( 31 ) 

k 3=1 3=1 

for all k = 1,2,..., N. 

We have seen before that, because of the U(l) gauge symmetry, each term in Vkim has an 
equal number of 0's and 0^'s, viz. the sum of all the columns of D vanishes identically. This 
means that D has a right-eigenvector (1, 1, 1) T corresponding to the eigenvalue 0. Now, 
D is a circulant matrix, because of the invariance under (151) . which means that each row of D 
is equal to the previous row moved one step to the right. As a consequence, the sum of the 
D-matrix elements over any column of D is equal for all the columns and it is equal to the sum 
of the D-matrix elements over any row of D, viz. to zero. Thus, the sum of all the rows of D 
vanishes identically. This means that D has a left-eigenvector (1, 1, . . . , 1) corresponding to 
the eigenvalue 0, or, equivalently, that D T has a right-eigenvector (1, 1, ... , 1) T corresponding 
to the eigenvalue 0. Therefore, the extremum points of J given by eqs. ( 131]) have 

sin -01 = sin ^2 = ■ • • = sin-^Tv- (32) 
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It must moreover be kept in mind that, since the sum of all the rows of D is identically zero, 

N 

ipj = mod 2tt. (33) 

j=i 

Equations ( |32l) and ( |33l) are the ones which yield the values of the ipj a t the minimum. 

Once the solutions of eqs. ( 1321) and ( 133]) for the are found, the phases 6k must be obtained 
from ijjj = djkdk- The matrix D has determinant zero and therefore we cannot invert this 
matrix equation. Rather we may, without loss of generality, set Q\ = and then define the 
(N — 1) x (TV — 1) matrix D which is identical to D with its first row and its first column both 
removed. We then have 

N 

^ = E( 5 U^ for j = 2,..., iV, (34) 

k=2 

which may be inverted: 

N 

0* = £P~%^ for k = 2, . . . ,N. (35) 
i=2 

Equation f )35|) is a non- homogeneous matrix equation for the 6^- Its general solution can be 
written as one particular solution plus the general solution of the corresponding homogeneous 
equation, which is 

N 

#M = E(^%( 27m ,)' (36) 
i=2 

with integer rij. Equation (13T)]) has exactly the same form as eq. and we already know that 
the set of its solutions is the rephasing symmetry group of the potential. 

When Afcim < the potential is minimized by maximizing J. This is simply done by setting 
all the ipj = 2imj with integer rij. Then eq. (135]) becomes identical with the homogeneous 
eq. (|36|) . This means that the 6k that we will get are essentially trivial, i.e. Ok = up to a 
transformation of the rephasing symmetry group G of Vk\ m - Hence, there is no spontaneous 
CP violation when \k\ m < 0. 

In the opposite case Xkim > 0, we need to minimize J instead. Now the analysis depends 
on whether iV is even or odd. In the case of even iV the relevant solution of eqs. (132]) and ( 133]) 
is ipj = 7i for all j = 1,2,..., N. This is a CP-conserving situation. Indeed, suppose that 
some non-trivial phases 6k solve eqs. (134"]) with all the ipj = ix. Then, the phases —6k also 
solve the same equations. Even if the two alignments 6k and — 6k do not coincide, they differ 
by a solution to the homogeneous eq. (I3T)]) . i.e. they are linked by a rephasing symmetry 
transformation. We then obtain that 

N 

«>;=EM*;>o> (37) 

3=1 

where U is a unitary symmetry of the potential. This means that the vacuum is invariant 
under a (generalized) CP transformation which is a symmetry of the Lagrangian. We have 
thus found that, in cases with a stable rephasing symmetry, there can be no geometrical CP 
breaking if N is even. 
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In the case of odd N, eqs. (1321) and (13"3"i) have solutions 



7T 



(38) 



which are CP-conjugate of each other and yield J = —Ncos(tt/N). Note that this result is 
insensitive to the exact form of the potential, namely to the values of k and m in Vkim- 

Given the solution (138]) for the ipj, we must construct the particular solution of eq. ( 135]) 
through 



A? 



e k = ^Y,i D ' l )k 3 forfc = 2,...,iV. 



(39) 



i=2 



This particular solution depends on the structure of the potential, viz. of the values of k, m, 
and N. However, often there exists another particular solution of eq. f )35|) . which has the 
simple form 

Ok = Vh = ^ {k ~ 1) for A; = 1, 2, . . . , N, (40) 

satisfying our Unsafe Q\ = 0. We stress that the existence of the solution (HOI) should always 
be explicitly verified in each particular case. 



3.3 First example: V213 for N = 3 
We first consider the case of the potential 



^phase — A213 
of ref. [7]. In this case 



(0102) (0103) + (0 2 3 ) (0 2 0l) + (030l) (030 2 ) + H.C.] (41) 



D 




and 



(42) 



(43) 



Therefore, according to eq. (j39j) . a particular spontaneously-CP-breaking solution with geo- 
metrical phases is ^2 = #3 = "00 while d\ = 0. The value of -0o is given by eq. ( )38|) : -00 = ±27r/3 
for the two CP-conjugate solutions. 

Notice that in this case the vector rj of eq. ( 140]) is not a solution to eqs. (|34|) . since 



— £> 

3 



does not coincide with ± (27r/3) (1, 1) as one would wish. 

The general solution involves the solution to eq. ( 1361) . which is identical with the solutions 
for the rephasing symmetry group. Since the potential (l4T|) is of the form V213 with N = 3, 





2tt 



(44) 
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the symmetry group is Zi S4 i = Z 3 with 2 — > ^02, 03 W 
CP-breaking solutions are 



h, d 2 , 9 3 ) 
h, O2, O3) 
h, #2) O3) 



0. 



2tt 2tt 

T' T 



47T 

0, T ,o 



0, 0, yj 



j. Therefore, the geometrically- 

(45a) 
(45b) 
(45c) 



and their CP conjugates. The solution (145aj) is equivalent, through a rephasing of the doublets, 
to 9 2 , 9s) = (— 2n/3, 0, 0) = (47r/3, 0, 0), which renders evident the symmetry among 
the three solutions. To the solutions (jlSjl must of course be added their CP conjugates. 



3.4 Second example: V213 for N = 5 

In this case the matrix D and its Smith normal form are 

-1 



/ 2 -1 -1 \ 

2-1-10 

2 -1 -1 

-1 2 -1 

\ -1 -1 2 / 



and 



/ 1 \ 

10 

10 

11 

\ ) 



(46) 



respectively. The symmetry group is therefore Zn = Zi S6 i which, according to section 12.21 is 
generated by phase rotations through 



^5 = ^ (0, 1, 10, 3, 6) . 



(47) 



To keep the notation short, we write here and below not the phase transformation themselves, 
but only their phase shifts; then, instead of taking powers of the transformations, we must 
write multiples of the phase rotations. The value of -00 is ipo = n ± 7r/5 = ±47r/5. Since 



D 



( 2 -1-1 \ 

2-1-1 

2 -1 

\ -1 2 / 



(48) 



the particular values of the 9k of the form (T40 



m = Y (0,1,2,3,4) 



(49) 



indeed provide a particular solution to eq. ( 134")) . since one obtains ip2 = 03 = "04 = 05 = 47r/5. 
Therefore, the generic calculable phases of the VEVs are in this case 



2tt 



9 (q) =rj s + qa B = — [(0, 11, 22, 33, 44) + q (0, 5, 50, 15, 30)] 



55 



(50) 
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for q = 0, 1, ... , 10. Including their CP conjugates, there are 22 distinct degenerate global 
minima with calculable phases, related among themselves by rephasing symmetry transforma- 
tions. 

Let us stress that any cyclic permutation of the phases in eq. ( |50|) also leads to viable 
minima. These are not additional minima, since they are already included in eq. ( 150]) . Indeed, 
we remind the reader that all the phases are defined up to overall phase rotations — this is 
the gauge symmetry U(l) discussed at the end of section [27T1 We always set Q\ = in order 
to pick a representative point of each family of relative phase configurations. Thus, a cyclic 
permutation of r] 5 is equivalent to r/ 5 . Moreover, a cyclic permutation of 05 may be expressed 
as a multiple of 05; for instance, if b is the generator of the cyclic permutation group Z 5 defined 
by the transformation ([5]), then 

h~ x a 5 b = — (6, 0, 1, 10, 3) = — (0, 5, 6, 4, 8) = 5a 5 . (51) 

Therefore, a cyclic permutation of any 0® m eq. (150]) just yields 6 (bq \ 

Let us also mention that, because of eq. ( ]5"T]) . the symmetry group of the potential is the 
semi-direct product Z n x Z 5 , plus the gauge-£7(l) symmetry group of overall phase rotations. 



3.5 Other potentials for N = 5 

In section 12.31 we have found four distinctly looking potentials, all of them protected by 
stable rephasing symmetries of the Z n type. They were generically written Vi + fc i i i i +2 fc for 
k = 1,2,3,4. In subsection 13.41 we have analyzed the case k — 1, i.e. the potential V213. 
Geometrical CP violation for the remaining three potentials can be easily written by using 
the same parameter k. Namely, for N = 5, the potential Vi+k, i,i+2fc has a particular solution 
pn^,, where p = 1 for k = 1 and k = 4 and p = 2 for k = 2 and k = 3. In each case, adding 
to this particular solution an arbitrary transformation from the corresponding rephasing sym- 
metry group generated by ( 12"3"]) yields the calculable phases of all possible global minima. We 
thus conclude that any potential protected with a stable symmetry for N = 5 can exhibit 
geometrical CP breaking, with the VEVs acquiring calculable phases. 

Noticing that 2r/5 is essentially equivalent to r/5 up to a permutation of the doublets, we 
conclude that in all cases the VEVs exhibiting geometrical CP violation for N = 5 represent 
the vertices of the regular pentagon on the complex planed 

3.6 Third example: V213 for N = 7 
The symmetry group is Z43 with generator 

a 7 = — (0, 1, 42, 3, 38, 11, 22) . (52) 

It can be verified again that 

777 = y (0, 1, 2, 3, 4, 5, 6) (53) 
is a particular solution, so that a generic global minimum has calculable phases 

= nr + qa 7 for q = 0, 1, . . . , 42. (54) 
3 The relevance of the regular pentagon to geometrical CP violation with N = 5 was also noted in ref. [9]. 
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For the other potentials protected by the stable symmetry, listed in ( 1211) as Vi + fc,i,i+2fe, we 
have the same situation as in the case N = 5, namely, in each case a particular solution is prj^ 
with p = 1 for k — 1, 6, p — 3 for k = 2,5, and p = 2 for = 3, 4. However, all these seemingly 
disparate cases are in fact, up to a permutation of the doublets, of the same form as 777. 



3.7 Fourth example: V214 for N = 7 
We now have 



D 



(55) 



/ 2 -1 -1 \ 

2-10-10 

2 -1 -1 

2 -1 

-1 2 -1 

\0-10 2/ 

-4, -4, -4, 3, 3, 10) T = (3, 3, 3, 3, 3, 3) T modulo 7. Thus, 



Therefore, D (1, 2, 3, 4, 5, 6) T = 
777 is once again a particular solution. The general solution is 

02, 03, 04, 05, 06, 7 ) = r?7 + ^ (0, 1, 8, 28, 23, 17, 4) 



(56) 



according to eq. (|25|) . Once again, there are further cases Vi-ffc, i, for k = 2, . . . , 6, but 
it turns out that they are all identical to the case V214 up to permutations of the doublets. 
Therefore, all the cases with VEVs exhibiting geometrical CP breaking for N = 7 represent 
the vertices of the regular heptagon on the complex plane. 



4 Beyond stable symmetries for even N 



In section 13.21 we have shown that in the case of an even number of doublets there is no room 
for geometrical CP violation when the potential is protected by a stable symmetry. One may 
ask whether this conclusion changes when the requirement of stability is removed. In the 
following we present an argument against this possibility. 
Let us firstly consider the simple case N 



4. As we have proved in section I2.3[ the 
potentials V213 and V234 are then invariant under the same rephasing group Z 5 . Suppose then 
that both terms are present in the potential: 



v; 



phase 



A 



+A' 



^3 



A<t>2 



+ cyclic + H.c. 



404 



+ 



403 



401 



+ H.c. 



(57) 



Using the Ansatz ((py = ue k , one obtains 



(Ve) = 2« 4 (AJ+A'J') 



(58) 



where 



J = COS Ipi + COS 1p2 + COS 1p3 + COS ^4, 
J' = COS (^1 + ^3) + COS (i> 2 + ^4) ■ 



(59a) 
(59b) 
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If only the first term with positive A was present, there would be no spontaneous CP violation, 
because the minimum would be given, up to a rephasing symmetry, by phases (0i, 2 , 3 , 4 ) = 
(0, 7T, 0, 7i"). The question is whether the presence of the second term may lead to non-trivial, 
calculable phases. 

To answer this question we have used the geometrical minimization methods described 
in ref. [TO]. Namely, we have considered the (J, J') plane, drawn the region covered by all 
possible phase configurations, and inspected it for the presence of sharp corners. One may 
find the shape of the allowed region analytically^ 

I J' I < 2, (60a) 
J 2 < 4J' + 8. (60b) 

This region has only two sharp vertices, (J, J') = (±4, 2). Both of them correspond to CP- 
conserving alignments: (6 U 6 2 , 63, 4 ) = (0, 0, 0, 0) for (J, J') = (4, 2) and (0 1; 2 , 3 , 4 ) = 
(0, 71, 0, 7i") for (J, J') = (—4, 2). All other alignments are not calculable. This proves the 
absence the geometrical CP violation in this model. 

The above arguments may be extended to a generic situation with even N. We start with 
some set of terms V^ m and define the corresponding J = J2j f =i cos ' l f J j- The fact that Vki m is 
not protected by a stable symmetry means that it is possible to write a second set of terms, 
which leads to a J' containing — at least for small N — the cosines of pairs of ip/s. One then 
derives the allowed region in the (J, J') plane: 

N 

\J'\ < -, (61a) 
N 2 

J 2 < NJ' + —, (61b) 

which has the same geometry as before and has only two CP-conserving corners. 

This discussion suggests the conclusion that for an even number of doublets it is much 
more difficult — maybe impossible — to achieve geometrical CP violation with renormalizable 
potentials. Then the only way to do it would be to include higher order (non-renormalizable) 
terms, as discussed in ref. [9]. 



5 Conclusions 

Since the fundamental origin of CP violation is still unknown, any theoretical framework in 
which CP violation arises naturally is worth a dedicated study. Geometrical CP violation, as 



4 Equation (|60b|) may be derived in the following way. Consider four planar vectors Hi, H2, H3, and H4, all 
of them with unit length and with polar angles ipi, ip2, —1P3, and —ipi, respectively. We define two vectors 
p = Hi + n,3 and q = + H4. Then, 

4 + 2J' = 2 + 2cos(<0i + ^ 3 ) + 2 + 2cos(<0 2 +^4) = |"i +n 3 | 2 + \n 2 + n 4 \ 2 = |pf + \qf . 

On the other hand, J is the projection of p + q on the horizontal axis. Then, 

J 2 < \P + qf = 2 (|pf + \qf ) - \p- qf < 2 (|pf + \qf ) = 8 + W . 
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suggested in ref. [TJ, is one such framework. The idea of calculable CP- violating phases which 
spontaneously arise in the vacuum of an extended Higgs sector is rather general and goes 
beyond the particular 3HDM example presented in ref. [TJ . In this paper we have investigated 
this phenomenon in the context of models with more than three Higgs doublets. We have 
shown how one may find the rephasing symmetry group of any cyclic iV-Higgs-doublet potential 
and have given several examples of geometrical CP violation for N = 5 and N = 7. In all 
those cases we have found a VEV alignment in the form of a regular polygon in the complex 
plane, distinct from what happens in the N = 3 model of ref. [TJ. In contrast, for an even 
number of doublets we have proved the impossibility of geometrical CP violation for phase- 
sensitive potentials which are protected by a 'stable' symmetry and have given an argument 
that suggests that it might be impossible to achieve geometrical CP violation even for more 
general renormalizable potentials. 
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